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We study the minimal supersymmetric standard model (MSSM) with minimal flavor violation 
(MFV), imposing constraints from flavor physics observables and MSSM Higgs searches, in light of 
the recent discovery of a 125 GeV Higgs boson by ATLAS and CMS. We analyze the electroweak 
vacuum stability conditions to further restrict the MSSM parameter space. In addition, a connection 
to ultraviolet physics is shown via an implementation of renormalization group running, which 
determines the TeV-scale spectrum from a small set of minimal supergravity parameters. Finally, 
we investigate the impact from dark matter direct detection searches. Our work highlights the 
complementarity of collider, flavor and dark matter probes in exploring the MSSM, and shows that 
even in a MFV framework, flavor observables constrain the MSSM parameter space well beyond the 
current reach of direct SUSY particle searches. 



I. INTRODUCTION 

Despite the null results for direct searches of supersym- 
metric particles at the LHC [1-29] , models of supersym- 
metry (SUSY) remain among the most well-motivated 
and popular extensions of the Standard Model (SM). Be- 
sides direct searches, there exist numerous ways to indi- 
rectly probe SUSY models, e.g. with low energy flavor 
observables, from dark matter direct detection results, 
and through Higgs properties. The discovery of a new 
particle at the LHC with a mass of ~ 125 GeV compati- 
ble with a SM-like Higgs boson [30, 31] has far reaching 
consequences for any model of New Physics (NP) with 
a non-standard Higgs sector. Indeed, the LHC Higgs 
results have motivated numerous studies of their impli- 
cations in the context of the Minimal Supersymmetric 
Standard Model (MSSM) and its variants [32-75]. 

A SM-like Higgs with a mass of M h ~ 125 GeV can 
be accommodated in the MSSM as long as stops are ei- 
ther very heavy or strongly mixed. Interestingly enough, 
large stop mixing unavoidably leads to irreducible con- 
tributions to low energy flavor observables, in particular 
in the Flavor Changing Neutral Current (FCNC) decays 
B s —> and B — > X s j, even if all soft masses are 

flavor blind. Correspondingly, rare B decays can be used 
to set non-trivial constraints on MSSM parameters. 

In this work, we discuss the status of the MSSM, 
in view of the recent Higgs search results from the 
LHC [30, 31] and the Tevatron [76], the recent strong lim- 
its on MSSM Higgs bosons in H/A -t bb searches [78, 79] 
and H/A — > t + t~ searches [77] , the latest results in B 
physics, in particular the recent evidence for B s — >• fi + n~ 
from LHCb [80], the updated results on B — >• tv from 
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Belle [81] and BaBar [82] as well as on B — > X s j from 
BaBar [83], and also the updated XenonlOO bounds on 
dark matter direct detection [84]. We will assume that 
the flavor structure of the SUSY particles is determined 
by the principle of minimal flavor violation (MFV) [85- 
88], i.e. we will assume that the SM Yukawa couplings 
are the only sources of flavor violation. This is moti- 
vated by the absence of any unambiguous deviation from 
SM expectations in flavor observables. We emphasize 
that even in this restrictive framework, flavor observ- 
ables play an important role in constraining the viable 
parameter space of the MSSM. Indeed, flavor bounds can 
be stronger than bounds from direct searches for SUSY 
particles in various regions of parameter space. This is 
particularly true for large values of tan (3, where loop- 
induced flavor changing couplings of the heavy Higgs 
bosons of the MSSM give enhanced contributions to 
FCNC processes. In the MSSM with large tan/3, di- 
rect searches of the heavy Higgs bosons also become es- 
pecially sensitive and, moreover, the exchange of heavy 
Higgs bosons can also lead to large dark matter direct de- 
tection cross sections, giving additional complementary 
means to probe this region of parameter space. 

In the large tan /3 regime, loop corrections to Higgs- 
fermion couplings can be significant and it is mandatory 
to resum tan /3-enhanced terms to obtain reliable predic- 
tions for any observables that depend on these couplings 
in the MSSM. We provide comprehensive analytical ex- 
pressions for all the relevant loop-corrected Higgs cou- 
plings, loop corrections to the SM-like Higgs mass, Higgs 
production and decay rates, contributions to flavor ob- 
servables, and dark matter direct detection cross sections, 
consistently taking into account the most general struc- 
ture of the soft SUSY breaking terms compatible with 
the MFV ansatz. In particular, we include effects from 
the bottom Yukawa coupling as well as the tau Yukawa 
coupling, as they are relevant for large tan/3. This goes 
beyond the analyses in [89-91], where bottom Yukawa 
effects in the squark masses were neglected. 
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We also give a detailed treatment of gaugino loop con- 
tributions to FCNC processes that can arise from a mass 
splitting between the left-handed squarks of the first two 
and the third generations. We highlight that in order 
to discuss the gaugino contributions to FCNC processes 
in the large tan/3 regime, both the squark mass split- 
ting and the alignment of this splitting with the quark 
Yukawas must be considered. 

Putting together all the relevant experimental con- 
straints coming from current Higgs, flavor and dark mat- 
ter sectors on the MSSM parameters, we point out re- 
gions of the MFV MSSM parameter space where these 
constraints are minimized. We also discuss the robust- 
ness of these bounds and to which extent they can be 
avoided. We take a phenomenological approach and 
treat the MSSM parameters as free parameters at the 
TeV scale: however, we augment this discussion with a 
study of renormalization group equation (RGE) effects, 
assuming minimal supergravity (mSUGRA)-like bound- 
ary conditions at a high scale and monitoring the generic 
spectrum of SUSY particles and their mass splittings in- 
duced by the running. 

In Sec. II, we review the MSSM with minimal flavor vi- 
olation in the quark sector. The impact of Higgs searches 
at the LHC and Tevatron on the MSSM are analyzed in 
Sec. III. We use the recent results indicating the pres- 
ence of a SM-like Higgs as well as dedicated searches for 
the additional Higgs bosons of the MSSM. In Sec. IV, we 
study constraints on large \i tan (3 from vacuum stability 
considerations. Constraints from B physics observables 
are analyzed in Sec. V. In Sec. VI, bounds on the MSSM 
from dark matter direct detection are considered. We 
conclude in Sec. VII. 



II. THE MSSM WITH MINIMAL FLAVOR 
VIOLATION 

In the following, we briefly review the MSSM with 
MFV. Throughout this work, in addition to MFV, we 
also assume minimal CP violation, i.e. the phase of the 
CKM matrix is the only source of CP violation, while all 
the MSSM parameters are CP conserving. We discuss 
the MFV structure of the sfermion masses in Sec. II A. 
Aspects of the Higgs spectrum that are relevant for our 
work are briefly reviewed in Sec. II B. In Sec. II C, we 
detail the tan /3-enhanced loop corrections to the Higgs- 
fermion couplings, allowing for the most general squark 
spectrum compatible with our MFV and CP conserva- 
tion assumptions. 



A. Sfermion Spectrum 

The soft SUSY breaking terms that give mass to the 
squarks, i.e. the soft masses, Wq, m 2 D and mfj, and tri- 
linear couplings, A d and A u , are possible sources of flavor 
violation. Generic flavor violating entries in these matri- 



ces are strongly constrained by flavor physics data. A 
simple approach to address this "SUSY flavor problem" 
is the principle of minimal flavor violation [85-88] , which 
states that the SM Yukawa couplings are the only sources 
of flavor violation even in extensions of the SM. In the 
context of the MSSM, this implies that the soft terms can 
be expanded in powers of the Yukawa couplings. In the 
super-Cabibbo-Kobayashi-Maskawa (super-CKM) basis, 
where squarks and quarks are simultaneously rotated 
to obtain diagonal Yukawa couplings, the soft masses 
are [85] 

m 2 Q = m 2 Q (l + b 1 V^y 2 l V + b 2 y 2 

+b 3 (y 2 d V^y 2 u V + V^y 2 u Vy 2 )) , 
™>u = ™u (l + b ivl) > 

m 2 D = m 2 D (l + b 5 y 2 ) , (1) 

where y u and yd are the diagonal up and down quark 
Yukawa matrices and V is the CKM matrix. The soft 
mass nig enters the left-left block of the down squark 
mass matrix, while VmqV^ enters the up squark mass 
matrix. The generic structure in (1) is always gener- 
ated by RGE running if flavor blind boundary conditions 
are assumed at a high scale [92, 93]. The parameters 
bi lead to splittings between the squark masses. To be 
specific, the parameters 64 and 65 generate a splitting 
between the masses of the first two and the third genera- 
tions of right-handed up and down squarks, respectively, 
while the parameters b\, b 2 , and 63 generate a splitting 
between the first two and the third generations of left- 
handed squarks. Note that the parameters b 2 , &3, and 65 
only become important for large values of tan /3, where 
the bottom Yukawa is 0(1). As we are particularly inter- 
ested in the large tan/3 regime, in the following we will 
take all the above masses as independent parameters and 
use rriQ 3 , m 2 D3 , and for the stop and sbottom masses 
and ttiq, m 2 ^, and vn\j for the masses of the first two gen- 
erations, which are degenerate to an excellent approxi- 
mation in this setup. This is analogous to the pMSSM 
framework [94] frequently studied in the literature. 

We stress, however, that the parameters b\, b 2 , and 
&3 also induce flavor violating entries in the left-handed 
squark mass matrices. These entries are proportional 
to small CKM angles and lead to controlled but non- 
negligible contributions to FCNC processes. In fact, due 
to SU(2)l invariance, the left-left blocks of the up and 
down squark mass matrices are related by a CKM rota- 
tion, and therefore any splitting in the diagonal entries of 
the left-handed soft masses rhg unavoidably induces off- 
diagonal entries in the up or down squark mass matrices. 
Moreover, distinct flavor phenomenology arises depend- 
ing on which of the b±, b 2 , or 63 parameters is respon- 
sible for the splitting. In particular, a splitting induced 
by 61 (62) is aligned in the up- (down-) sector and will 
only lead to off-diagonal entries in the down (up) squark 
mass matrix. The parameter 63 induces flavor violation 
in both the up and down squark masses matrices. All 
flavor observables that we will discuss in the following 
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depend on the combination b\ 
one additional parameter 



bay 2 . We thus introduce 



c 



hy 2 + b 2 y 2 b + 2b 3 y 2 y 2 



(2) 



which reflects the alignment of the splitting in the left- 
handed squark masses and hence parametrizes the frac- 
tion of the splitting in the masses leading to flavor vi- 
olation in the down sector. We assume C is real in the 
following. 1 We see that formally ( — 1 + 0(y 2 ). If we 
consider a splitting in the squark masses that is radia- 
tively induced through RGE running, then considering 
only the top Yukawa in the running leads to £ = 1. Bot- 
tom Yukawa effects become important for large tan j3 and 
can lead to < C < 1- Typically we expect that y b is at 
most as large as yt, however, which implies 1/2 < £ < I. 

We note that an expansion analogous to (1) also ex- 
ists for the trilinear couplings [85]. In particular, higher 
order terms in the expansion can lead to flavor violating 
trilinear terms. Such terms only lead to corrections of the 
holomorphic Higgs couplings, however. These corrections 
can induce flavor changing neutral Higgs couplings, that 
are especially interesting beyond MFV, where the corre- 
sponding effects can be chirally enhanced [95, 96]. In the 
MFV framework considered here, these effects are less 
important compared to contributions that are related to 
the loop-induced non-holomorphic Higgs couplings. The 
only relevant trilinear couplings for our analysis are those 
for the third generation squarks, A t and At,, which we will 
take to be independent parameters. 

For simplicity, we will also assume universal soft 
masses m 2 L and m 2 E , in the slepton sector. The phe- 
nomenology of flavor non-universalities in the lepton sec- 
tor will be reserved for future study. The only relevant 
trilinear term in the slepton sector is the tau trilinear 
coupling A Tl which, along with A t and A b and all other 
parameters, we will take to be real. 



B. Higgs Spectrum 

The physical Higgs spectrum of the MSSM consists of 
two neutral scalar bosons h and H, one neutral pseu- 
doscalar A, and a pair of charged Higgs bosons H ± . At 
tree level, the full spectrum is determined by only two 
real parameters: the mass of the pseudoscalar Higgs, 
Ma, and the ratio of the two vacuum expectation val- 
ues, tan/3 = tp = v u /v d , with v 2 +v^ — v 2 = 17 '4 2 GeV 2 . 
In the so-called decoupling limit, M\ ^> M^, the masses 



of the Higgs bosons, A, H and H , are 



Ml 



M\ 



M 2 H± 



M\ 



(3) 



1 Note that while b\ and 62 have to be real due to hermiticity of 
the squark masses, 63 can in principle be complex. Indeed, as 
shown in [93] , a tiny phase for ^3 is always generated during RGE 
running. 



In this limit, the mass of the lightest Higgs h is given at 
tree level by 



Mi ~ M 2 Z cos 2/3 



(4) 



As is well known, moderate or large values of tan /? and 
large 1-loop corrections are required to lift Mh up to phe- 
nomenologically viable values. Moreover, at large tan /3, 
the sbottom and stau 1-loop corrections can lower M% by 
a few GeV, which cannot be neglected given the current 
Higgs mass precision data. The dominant stop, sbottom, 
and stau loop contributions for large tan/3 read 
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where X t = A t — /i/tan/3 w A t for large tan/3, and 
m^, m- h and nif are the average stop, sbottom, and stau 
masses, respectively. The stop loop corrections, reported 
in the first line of (5), are maximized for A t ~ \f&mj.. 
The contributions from the sbottom and stau loops, in 
the second and third lines, always reduce the light Higgs 
mass and can be particularly important for large tan /3, 
large values of the Higgsino mass parameter, /1, and light 
sbottom or stau masses [97]. The factors come from 
an all-order resummation of tan/3 enhanced corrections 
to the Higgs-fermion couplings and are discussed in detail 
in Sec. II C. 

The couplings of the lightest Higgs to SM fermions 
and gauge bosons are mainly controlled by tan /3 and the 
angle a that diagonalizes the mass matrix of the two 
scalar Higgs bosons. If 

a = /3-7r/2, (6) 

the couplings of h are exactly SM-like. At the tree level, 
Eq. (6) holds up to corrections of order M^/{tpMl). 
Correspondingly, for large tan/3 and moderately heavy 
Ma, the couplings of h are already SM-like to a good 
approximation. At 1-loop, Eq. (6) gets corrected by an 
additional term ~ A7V 2 jM 2 A , where A7 is a loop- induced 
quartic Higgs coupling that reads 

3 m\\iA 

A 7 — 7^ 9 A 9 
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If A7 is sizable, corrections to the light Higgs couplings 
can become relevant, as discussed in the next section, and 
are constrained by the SM Higgs searches at the LHC. 
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C. Higgs Couplings to Fermions 

At tree level, the MSSM Higgs sector is a 2 Higgs dou- 
blet model of type II, where only H u couples to right- 
handed up quarks and only Hd couples to right-handed 
down quarks and leptons. The Yukawa interactions thus 
have the following form 



•^Yuk — (Uu)ij 



{Vcdij H.iQ.l), 

+ h.c. . 



(8) 



As a consequence, the couplings of the neutral Higgs 
bosons to fermions are flavor diagonal in the mass 
eigenstate basis. At the loop level, however, "wrong" 
Higgs couplings are generated and lead to potentially 
large threshold corrections to the masses of down type 
quarks and leptons [98-101] as well as CKM matrix el- 
ements [102]. They also significantly modify charged 
Higgs couplings to quarks [103, 104] and generate fla- 
vor changing neutral Higgs couplings [95, 96, 105-110]. 
All these effects become particularly relevant in the large 
tan/3 regime, where the inherent 1-loop suppression can 
be partly compensated. In the following, we analyze 
the form of the neutral and charged Higgs couplings 
with fermions in the phenomenologically motivated limit, 
v 2 <C M£ usy (see [110] for a discussion of the regime 
v 2 ~ A^susy)- We consistently take into account the 
most generic MFV structure of the squark masses as dis- 
cussed in Sec. II A. In particular, we consider splittings 
between the first two and the third generation squarks in 
the left-handed as well as right-handed sectors. 

Once the 1-loop corrections are taken into account and 
we have diagonalized the quark mass matrices, the cou- 
plings of the neutral Higgs mass eigenstates to fermions 
have the following generic form 

L ^ D E i feV + $< H + %M ( 9 ) 

9,9' 

+ J2 {$h + $H + itf A) + h.c. , 



where we have neglected flavor changing couplings to lep- 
tons, which are not relevant for our analysis. Using the 
notation £ qq = and again for large tan/3, the flavor 
conserving couplings of the heavy scalar and pseudoscalar 
Higgses, normalized to their respective SM Yukawas, are 



K q -K q - ^ e q 



for q = u,c,t , (10) 



iovq = d,s,b, (11) 

q q 1 + e a t B 

(12) 



pH ^ pA 
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In the above expressions, tan /^-enhanced corrections to 
the couplings are resummed to all orders and the fac- 
tors parametrize the loop induced "wrong" Higgs cou- 
plings. The exact form of each e t in terms of MSSM 
parameters is discussed at the end of this section. Since 



we assume the MSSM parameters to be CP conserving, 
all e, parameters are real. 

Among the flavor changing Higgs couplings only the 
coupling of a right handed bottom with a left-handed 
strange quark will be relevant in the following discussion. 
Normalized to the bottom Yukawa of the SM, we have 
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(l + e 6 t /3 )(l + eoi /3 ; 



V tb V t 



tbVts i 



(13) 



where eo is defined as eo = e& — £fcj an d £fc is discussed 
in detail below. 

The couplings of the light Higgs boson, h, are exactly 
SM-like in the decoupling limit: = f£ = 1 and (* fq = 
for q' 7^ q. While non-standard effects in the couplings to 
up-type quarks are generically tiny even away from the 
decoupling limit, corrections to the couplings with down- 
type quarks and leptons decouple very slowly and can be 
relevant. We have 



Sq 1 ~ €f/t a 
Cp 1 + tftp 



for / = d, s, b, 



(14) 



The couplings of the physical charged Higgs bosons to 
fermions can be written as 



9,9' 

e 

For the couplings relevant to our analysis, we have 



(15) 



V tb 1 + e b t/3 
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1 + eetp 
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(16) 
(17) 
(18) 



where Vy are CKM matrix elements and e — e t — e FC . 
The parameter e' FC is the up-sector analogue of epc- 

As already mentioned, the various e factors in 
the above expressions parametrize loop-induced non- 
holomorphic Higgs couplings. They arise from Higgsino- 
squark loops, gluino-squark loops and wino-sfermion 
loops. We do not explicitly state the typically negligible 
contributions coming from bino-sfermion loops; however, 
they are included in our numerical analysis. 

For the bottom quark, we can decompose e b = e? + 



, where these contributions are 
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\xA t g([i 2 \m 2 Q ^m 2 V3 ) 



(19) 
(20) 
(21) 
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The loop function g is listed in the appendix, and has 
units of (GeV)~ 2 . Hence, the e factors generally exhibit 
non-decoupling as the SUSY mass scale increases. In 
particular, rescaling all the SUSY mass parameters, i.e. 
the squark masses, gaugino masses, the Higgsino mass 
parameter and the trilinear coupling by a common, ar- 
bitrarily large factor, leaves the e parameters invariant. 
For a degenerate SUSY spectrum with mass rh, we ob- 
tain g(m 2 ,m 2 ,m 2 ) = l/2m 2 . Our sign convention is 
that the left-right mixing entries in the top and bottom 
squark mass matrices are given by m t (A t — fj, cot 0) and 
m b {A b — /^tan/3), respectively. Furthermore, the gluino 
mass M3 is always positive in our convention. 

For the strange and down quarks, the Higgsino contri- 
bution is highly suppressed by small Yukawa couplings 
or CKM angles and only the gluino and wino loops are 
relevant: e S: d = ef d + e^d, where the e* d can be easily 
obtained from the corresponding e\ expressions, replac- 
ing third generation squark masses with second or first 
generation squark masses. 

For leptons, only the wino (and the bino) loops give 
contributions, and ef is given by ef with the sbottom 
masses replaced by the slepton masses. 

In case of the top quark, analogous to the bottom 
quark, we consider the gluino, wino, and Higgsino con- 



tributions: e t = ef + ej + ej . The expressions for ef 
and ef are trivially obtained from the corresponding e l b 
by replacing the relevant squark masses. The Higgsino 
contribution is explicitly given by 



t 2 



H _ _n\ 

* " 4ir2M* v {l + e b t p ) 2 



/j,A b g(v 2 ,m 2 Q 3 ,Tn 2 D3 ) 



(22) 

Here, ef is suppressed by the bottom quark mass and 
only becomes relevant for large values of tan /3. 

The flavor changing couplings, epc and e FC , can be 
decomposed as 



£FC = ( b + C e FC + C e 



C FC ' 



e' FC = ef + (1-C)e^ + (1-C)e 



w 

FC > 



(23) 
(24) 



with 



9 as 8 JLJT 



x (g(M 2 , m 2 Q3 ,m 2 D3 ) - g{Mlm 2 Q ,m 2 D3 )) ,(25) 
x (g(v 2 ,M 2 , m 2 Q3 ) - g(v 2 ,Ml,m 2 Q )) , (26) 



W "2 3 



and e'p C is obtained from e FC by replacing the right- 
handed sbottom mass, m D3 , with the right-handed stop 
mass, mu 3 . The ef and ef expressions were already 
given above. Note that m]j 3 enters both loop functions 
in (25) and hence, in general, e FC ^ e 9 b — ef , in contrast 
to the case where all right-handed down squarks have the 
same mass, mr> 3 — mo- Clearly, a splitting between the 



third and the first two generations of left-handed squarks 
induces non-zero e FC , e^ c and/or e'p C . 

As already described in Sec. II A, £ parametrizes the 
alignment of the left-handed squark mass with the quark 
masses. The case £ = 1 corresponds to a tuq that is 
aligned in the up sector such that the mass splitting be- 
tween the first two and the third generations leads to 
off-diagonal entries only in the down squark mass ma- 
trix. This in turn implies maximal flavor changing gaug- 
ino loop corrections to the Higgs-down quark couplings. 
The case C = corresponds to alignment in the down sec- 
tor, with no off-diagonal entries appearing in the down 
squark mass matrix. Generically, from RGE running, we 
expect 1/2 < C < 1. 



III. HIGGS COLLIDER SEARCHES 
A. SM-like 125 GeV Higgs 

The LHC experiments, ATLAS and CMS, have re- 
cently discovered a new particle with a mass of about 
125 GeV [30, 31]. This discovery is based on results from 
SM Higgs searches in the 77, ZZ and WW channels. 
The observed signals indicate that the new particle is a 
boson with spin or 2, and overall, they are in reason- 
able agreement with expectations for a SM Higgs. Other 
searches in the t + t - and bb decay channels are also be- 
ing pursued, but more statistics are needed in order to 
make conclusive statements. 

The most visible feature of the extracted signal 
strength in all the different channels under study is an en- 
hancement in the 77 decay rate in comparison to the SM 
rate. The decay rates into WW and ZZ gauge bosons 
are consistent with the SM values at the la level. The 
present experimental uncertainties in the signal strength 
in the various production and decay channels allow for 
many new physics alternatives. In particular, within su- 
persymmetric extensions, it is possible to enhance or sup- 
press the gluon fusion production with light stops, de- 
pending on the amount of mixing in the stop sector. It is 
also possible to suppress gluon fusion with light sbottoms 
that have sizable mixing driven by large values of /i tan /3. 
In all cases, enhancement of gluon fusion implies a sup- 
pression of the h — > 77 decay rate, and vice-versa. The 
overall effective gg — > h — > 77 rate, however, is governed 
by the enhancement or suppression of the gluon fusion 
production cross section. 

To achieve a net enhancement in the h — > 77 rate, 
uncorrelated with a simultaneous enhancement in the 
h WW/ZZ rates coming from an enhanced gluon fu- 
sion production or reduced h — > bb partial width, the exis- 
tence of new, light, charged colorless particles running in 
the loop is required. In the MSSM, the only two options 
are charginos, which only contribute for tan /? ~ 1 (dis- 
favored by a 125 GeV Higgs mass), and light staus with 
large mixing, i.e. large /ztan/3. A detailed discussion 
of the possible deviations from SM values of the produc- 
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tion and decay rates for a SM-like Higgs in the MSSM 
can be found in [39, 58]. Possible correlations with flavor 
observables have very recently been studied in [75]. 

While it is very interesting to investigate deviations 
from SM expectations in Higgs data that would point to- 
wards new SUSY particles within the reach of the LHC, 
we take a different approach in this work by assuming 
a Higgs boson with approximately SM-like properties. 
We concentrate on possible signatures of new physics 
that may appear in B physics observables, direct non-SM 
Higgs searches and dark matter direct detection searches 
within the MSSM with MFV, while fulfilling the require- 
ment of a 125 GeV SM-like Higgs. In this way, we show 
indirect effects from SUSY particles in flavor and Higgs 
physics in regions of parameter space beyond the present 
reach of the LHC. 



B. Searches for Heavy Scalars and Pseudoscalars 

Searches for the heavy neutral Higgs bosons of the 
MSSM have been performed in the H/A — > bb and 
H/A -> t+t~ channels both at the Tevatron [111-114] 
and the LHC [77-79, 115, 116]. 

Searches also exist for light charged Higgs bosons in 
top decays at both the Tevatron [117, 118] and the 
LHC [119-121]. For the MSSM scenarios considered in 
this work, however, the corresponding bounds are not 
competitive with the bounds from searches of the neu- 
tral Higgs bosons. 

In the large tan/3 regime, the cross sections for the 
heavy scalar and pseudoscalar Higgses rescale according 
to 



a bb^h x 



tt. SM w , 

<J a X 

gg->h \ £ 



(l + £btp) 2 
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tb, SM 1 — e tt/3 



bb, SM 



(27) 



,(28) 



evaluated at a common mass for all Higgs bosons. For 
large tan/3, the (Jbb-yH/A production cross sections can 
dominate over gluon fusion. We use HIGLU [122] and 
bbh@nnlo [123] to compute the respective SM cross sec- 
tions cr^™ and af^ h at the LHC. 

The most important decay modes of the heavy Higgs 
bosons are H,A — > bb and H,A — > t + t~. The corre- 
sponding partial widths can be written as 



Fflbb — ^Abb 
^Htt — L 'att 



t 2 

r SM H 
hbb (1 + ebttf 



hrr 



f 2 



(29) 
(30) 



where r|%, are the corresponding decay widths of a Higgs 
boson with the same mass as H and A and with SM-like 
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FIG. 1. Constraints in the M^-tan/J plane from direct 
searches of the neutral MSSM Higgs bosons at CMS and AT- 
LAS. The solid, dotted and dashed lines correspond to sce- 
narios (a), (b), and (c) as defined in Tab. I. The blue (green) 
regions are excluded by searches in the t + t~ (bb) channel. 



Scenario 


(a) (b) (c) (d) (e) 


H [TeV] 
sign (At) 


1 4 -1.5 1 -1.5 

+ + + - - 



TABLE I. Illustrative MSSM scenarios discussed in the text. 
All sfermion masses are set to a common value 2 TeV, the 
gaugino masses to 6A/i = 3M2 = M3 = 1.5 TeV. The trilinear 
couplings A t = At = A r are set such that the lightest Higgs 
mass is Mh = 125 GeV. 



couplings to bb and t + t . In our numerical analysis, we 
compute rlf f using HDECAY [124]. 

Note that the main dependence of the production cross 
sections and branching ratios is on tan/3 and the heavy 
Higgs masses. Dependence on other MSSM parameters 
enters only at the loop level through the tan /3 resumma- 
tion factors Cj. 

In our framework, the most important constraints 
come from the CMS bounds in the t + t~ channel [77], 
which are available up to masses of Ma = 800 GeV 
and the bb channel [78, 79] which cover heavy Higgs 
masses up to Ma < 350 GeV. Our estimates for the 
excluded regions from the H/A — > bb searches are shown 
in Fig. 1 in yellow-green and labeled with bb. We set 
all sfermion masses to 2 TeV and the gaugino masses to 
6M1 = 3M 2 = M 3 = 1.5 TeV. The solid, dotted and 
dashed contours correspond to a Higgsino mass parame- 
ter fi = 1 TeV (scenario a), 4 TeV (scenario b) and —1.5 
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TeV (scenario c), respectively. For every point in the 
M^-tan /3 plane, the trilinear couplings At = Ab = A T 
are positive and chosen such that the lightest Higgs mass, 
computed using FeynHiggs [125], is Mh — 125 GeV. 2 
The respective choices for fi lead to representative values 
for the tan/3 resummation factors of c^+e^ — 3.3x 10~ 3 , 
1.5 x 1(T 2 and -5.1 x 1CT 3 . 

As is well known, the bounds in the Ai^-tan /3 plane 
coming from the t + t~ channel are robust against varia- 
tions of the MSSM parameters. Indeed, the dependence 
of the production cross section on e& is largely cancelled 
by the corresponding dependence of the BR(^4, H —> 
t + t~) [126, 127]. In Fig. 1, we therefore simply re- 
port in blue the tan/3 bounds obtained in [77] in the 
so-called M™ ax scenario. We checked explicitly that the 
constraints are largely independent of the scenarios in 
Tab. I. We find that the constraints can only be weak- 
ened mildly for large tan j3 and Ma if the MSSM param- 
eters are such that tb is sizable and positive, as in sce- 
nario (b) 3 . We note however, that in the region with low 
tan /3, the bounds do depend to some extent on the SUSY 
spectrum, in particular the neutralino and chargino spec- 
trum. Indeed, for low tan/3, the heavy scalar and pseu- 
doscalar Higgs bosons can have sizable branching ratios 
in neutralinos or charginos if these decays are kinemati- 
cally allowed. The M™ ax scenario considered in [77] con- 
tains light neutralinos with M Xl ~ 95 GeV. For small 
tan /3, the obtained bounds from the searches in the t + t~ 
channel are therefore slightly weaker compared to scenar- 
ios with heavier neutralinos. 

The CMS searches in the bb channel [78, 79] are not yet 
competitive with the r + r~ searches, but might become 
important for large Ma in the future [128]. Compared 
to the t + t - searches, the bounds coming from the bb 
searches show a stronger dependence on the remaining 
MSSM parameters [126, 127]. In particular, for large neg- 
ative /i, the bounds become significantly stronger, while 
for large positive /Lt, the bounds can be weakened consid- 
erably. Note that for large negative [i and large tan/3, 
however, constraints from vacuum stability and pertur- 
bativity of the bottom Yukawa have to be taken into 
account. 

Since the theoretical precision of the light Higgs mass 
prediction in the MSSM allows shifts of a few GeV, 
we checked the extent to which the H/A — > t + t~ and 
H/A — > bb bounds depend on the exact value of the Higgs 



2 The Higgs mass, M^, is not a monotonic function in At and for 
a given sign of At there are typically two choices of At that lead 
to Mh = 125 GeV. We always take the At that is smaller in 
magnitude. 

3 The CMS results in the t + t~ channel are only available as con- 
straints in the M^-tan f) plane for the M ™ ax scenario. We trans- 
late these constraints into bounds on the corresponding trxBR 
and then reinterpret the cross section bounds as constraints in 
the M^— tan /3 plane for various choices of the other MSSM pa- 
rameters summarized in Tab. I. We assume constant efficiencies 
throughout this procedure. 



mass assumed in our analysis, Mh = 125 GeV. We find 
that varying the light Higgs mass in the range 122 GeV 
< Mh < 128 GeV does not change the constraints from 
H and A searches significantly. 



IV. VACUUM STABILITY 

Independent of experimental searches, large values of 
\x can be constrained based on vacuum stability consid- 
erations, particularly if tan (3 is also large. Indeed, large 
values of /itan/3 can lead to charge and color breaking 
minima in the scalar potential of the MSSM [129]. 

The trilinear couplings of the up- type Higgs, H u , with 
sbottoms, &l and Sr, and staus, ?l and ?n, are controlled 
by 

m T ^tan/3 
CD Vl + er tan/3 (i7 « TiTfl) + 



/Ztan^ f rjQT*7 s . , 

-, n\H u b L b R ) + h.c. 

v l + £fctan/3 v u L ' 



(31) 



For trilinear couplings that are large compared to the 
sbottom or stau masses, minima with non-zero vevs for 
the sbottom and/or stau fields can arise in addition to 
the standard electroweak minimum. If these minima are 
deeper than the electroweak minimum, the electroweak 
minimum becomes unstable and can decay. The corre- 
sponding regions of parameter space are only viable as 
long as the lifetime of the electroweak minimum is longer 
than the age of the universe. This corresponds to requir- 
ing that the bounce action, B, of the tunneling process 
is B > 400 [130, 131]. 

Our phenomenological flavor analysis is largely inde- 
pendent of the values of the stau masses, and non-zero 
stau vevs can always be avoided if rrif L and m? H are large 
enough. Nonetheless, we will consider a scenario where 
the squark and slepton masses are the same order and 
thus include both sbottoms and staus in the following 
analysis. 

Starting with the MSSM scalar potential, we restrict 
ourselves to terms that contain only the up-type Higgs, 
sbottoms, and staus, which are the degrees of freedom 
most relevant for large //tan/3. We consider three cases: 
(i) only terms with the up-type Higgs and staus, (ii) only 
up-type Higgs and sbottoms, and (iii) up-type Higgs, 
staus and sbottoms simultaneously. In each case, we 
search for additional minima in field space and estimate 
the bounce action for tunneling from the electroweak 
minimum into the deepest minimum of the potential. In 
the end, we apply the strongest of the three bounds. 

For each case, the second vacuum generally has sepa- 
rately nearly degenerate stau vevs and nearly equal sbot- 
tom vevs. In the case of the sbottoms, this is expected 
from the SU (3) D-terms in the scalar potential: 



C d 



\h\ 2 



\b R \ 2 



(32) 
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FIG. 2. Constraints from vacuum stability in the /i-tan /3 plane. We set the sbottom and stau soft masses to 2 TeV and the 
gaugino masses to 6Mi = 3M2 = M3 = 1.5 TeV. In the left (right) plot, the trilinear coupling of the stops is A t — 2 TeV 
(At = —2 TeV). The labeled contours show the values of the bottom Yukawa coupling. In the light red (light gray) regions, 
a charge and color breaking vacuum exists that is deeper than the electroweak breaking vacuum, but the electroweak vacuum 
has a lifetime that is longer than the age of the universe. In the dark red (gray) regions, the electroweak vacuum is not stable 
on cosmological time scales. Finally, in the black regions, one of the sbottoms becomes tachyonic. 



We can clearly see that at least for the 3-dimensional pa- 
rameter space in case (ii), deviations from equality of the 
squark/slepton fields along the path chosen to compute 
the action would come at the expense of large contribu- 
tions from the D terms. Therefore, to obtain an analyti- 
cal estimate for the bounce action, we consider a straight 
path in field space connecting the electroweak minimum 
and the charge and/or color breaking minimum. 

We then approximate the potential along the straight 
line by a triangle and use the analytical expressions 
in [132] to calculate the bounce action. We con- 
struct the triangle such that for a few chosen parame- 
ter points, the obtained bounce action agrees approxi- 
mately with the bounce action from the analytic expres- 
sion of the potential solved numerically by a standard 
overshoot /undershoot method. We further crosschecked 
our results with CosmoTransitions [133] taking into ac- 
count the up-type Higgs, the down-type Higgs, the sbot- 
toms, and the staus. Overall, we find good agreement 
with our approximate analytical approach. 

The constraints thus derived in the /Lt-tan/3 plane are 
shown in Fig. 2. We fix the SUSY masses as in the 
scenarios considered above, namely we assume degen- 
erate sfermion masses with rh = 2 TeV and gaugino 
masses with 6 Mi = 3M 2 = M 3 = 1.5 TeV. The tri- 
linear couplings we set to A t = 2 TeV in the left and 
to A t = —2 TeV in the right plot. In the white region, 
the electroweak minimum is the deepest minimum in the 
potential and therefore absolutely stable. In the light red 
(light gray) region, a charge (and possibly color) break- 



ing minimum exists that is deeper than the electroweak 
minimum, but the electroweak minimum has a lifetime 
longer than the age of the universe. In the dark red (gray) 
region the lifetime of the electroweak minimum is shorter 
than the age of the universe. Finally, in the black region, 
one of the sbottoms is tachyonic. 

The solid lines labeled in the plots show contours 
of constant bottom Yukawa couplings in the /i-tan /3 
plane. For large and negative /itan/3, close to the region 
where one of the sbottoms becomes tachyonic, the bot- 
tom Yukawa coupling becomes non perturbatively large. 

We observe that large negative values for (i are strongly 
constrained by the requirement of vacuum stability. This 
is because the tan/3 resummation factor, e^, in (31) is lin- 
early proportional to \i. It increases the trilinear coupling 
of the up-type Higgs with sbottoms for negative values 
of /i and can lead to a deep second minimum mainly in 
the field direction of the sbottoms. In particular, we find 
that values of p, tan /3 negative and large enough that the 
bottom Yukawa changes its sign (the parameter space 
in the upper left corner of the plot beyond the region 
excluded by tachyonic sbottoms) are excluded by the re- 
quirement of vacuum stability. For positive values of /1, 
the coupling of the up-type Higgs with sbottoms is re- 
duced while its coupling with staus is slightly enhanced 
by the e T term. In this region of parameter space, con- 
straints come typically from a second minimum in the 
stau direction. Positive values for /1 are less constrained 
than negative ones, and the allowed region for /i can be 
extended above /1 > 10 TeV for sufficiently heavy staus. 
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The viable regions of parameter space can be enlarged 
slightly when we allow for a splitting between the masses 
of the left- and right-handed sbottoms and staus. Nev- 
ertheless, we do not find any regions of parameter space 
where both vacuum stability and e b tan f3 < — 1 (which 
flips the sign of the tree level bottom Yukawa and hence 
changes the typical sign of the SUSY contribution to B 
observables) can be achieved simultaneously. In the end, 
we see that the scenarios discussed in the previous section 
are all compatible with bounds from vacuum stability. 



V. B PHYSICS OBSERVABLES 

Flavor observables play a crucial role in determining 
viable regions of parameter space of SUSY models. This 
is true both under the MFV assumption [89-91, 134-140] 
and if new sources of flavor violation are allowed [141- 
150]. 

Of particular importance in the MFV setup are rare 
B decays that are helicity suppressed in the SM, because 
SUSY contributions to these decays can be enhanced by 
tan/3 factors. Interesting processes include the tree level 
decay B — > tv, the purely leptonic decay B s —> 
and the radiative decay B — > X s j. Additional con- 
straints on the SUSY parameter space can be also de- 
rived from the (g — 2) of the muon. The (g — 2) M bound 
becomes particularly important if sleptons are only mod- 
erately heavy, which is a scenario that we do not consider 
here. 



A. B — y tv, B — y D^tv and K ->• /Ltf 

The decay B — > tv is a sensitive probe of extended 
Higgs sectors as it can be modified by charged Higgs 
exchanges at tree level [151]. The most important in- 
puts for the SM prediction are the CKM element \V u b\ 
and the B meson decay constant. Using the PDG value 
| Kb | = (3.89 ± 0.44) x 10" 3 [152], a conservative av- 
erage over direct determinations from inclusive and ex- 
clusive semi-leptonic B decays, and an average of re- 
cent precise lattice determinations of the decay constant 
f B = (190 ± 4) MeV [153-156], we find 



BR(B -> tv)sm = (0.97 ± 0.22) x 10" 



(33) 



While previous experimental data gave values for the 
branching ratio more than 2a above the SM prediction, 
a recent result from Belle [81] has a much lower cen- 
tral value. An average of all the available data from 
BaBar [82, 157] and Belle [81, 158] gives 



BR(B -> tv) cxp = (1.16 ± 0.22) x 10~ 4 



(34) 



This value is in very good agreement with the SM but 
still leaves room for NP contributions. 

Closely related decay modes that are also sensitive to 
charged Higgs effects are the B — > Dtv and B — > D*tv 



decays [159-163]. While predictions of the corresponding 
branching ratios suffer from large hadronic uncertainties 
coming from the B — » D and B — > D* form factors, 
the ratios BR(B -> Dtv)/BR(B -> Div) and BR{B -> 
D*tv)/BR(B — > D*£v), where £ = e or fj,, can be pre- 
dicted with reasonable accuracy in the SM [162, 164]. 
Interestingly, recent results from BaBar [165] on these 
ratios are around 2a above the SM predictions in both 
decay modes. Older results from Belle [166] give similar 
central values but with much larger uncertainties. 

Another interesting observable in this context is 
R^23 [167] that probes the tree level charged Higgs ex- 
change in the K — > jiv decay. The much smaller sensitiv- 
ity of K — > [iv to charged Higgs effects compared to the 
B decays is compensated by its extremely high experi- 
mental precision and the excellent control on theoretical 
uncertainties giving [167] 



R 



fj.23 



0.999 ±0.007 



(35) 



All the mentioned tree level decays depend in similar 
ways on possible new physics contributions in the MSSM 
with MFV. Defining 



X 



1 



B(K) 



f2 



M 2 H± (1 



£o(s)i/3)(l + Wtp) 



(36) 



we can write 



Rbtv 



(37) 



Rdtv — 



Rd-tv — 



Rfj.23 — 



BR{B -> tv) 
3R(B ->• tv) sm 

l-m 2 B+ X 2 B y , 

BR(B -> Dtv) 
BR(B Dtv) sm 

1-1.5 m T m b X 2 B + 1.0 m 2 T m 2 b X% ) , (38) 

BR{B -> D*tv) 



BR{B -> D*tv) sm 

(l-0.12m T m b X 2 B +0.05m 2 T m 2 b X B ) ,(39) 

BR(K -> \iv) 
BR(K -> hv)sm 

(l - m 2 K+ X 2 K 



(40) 



In Fig. 3 we show these ratios as function of Xb,k = 
X B K \ both for positive X 2 (solid lines) and nega- 
tive X 2 (dotted lines) in comparison with the experi- 
mental la and 2a ranges (dashed bands) from (33)-(35) 
and [162]- [165]. Here, positive X 2 illustrates destructive 
interference with the SM, while negative Xf illustrates 
constructive interference with the SM. 

We observe that agreement of theory and experiment 
in all three B observables is impossible to achieve. In 
particular the tensions in B — > Dtv and B — > D*tv 
cannot be addressed in the context of the MSSM with 
MFV, but require more radical approaches [168-172]. 
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FIG. 3. Branching ratios of the decays B — > TV, B — > Dtv, 
B — ¥ D*tv, and K — > fiv as functions of Xb, Xk = \J\X% K \ 
as appropriate and which are defined in the main text. The 
dashed bands show the la and 2a experimental ranges. The 
solid (dotted) lines are the theory predictions for positive 
(negative) X 2 giving destructive (constructive) interference 
with the SM amplitude. 

Considering MSSM contributions to K — > /iv and 
B — > tv, we observe that generally, Xb and Xk are equal 
to a good approximation. The only way to induce a dif- 
ference is through a splitting between the right-handed 
strange squark mass and the right-handed bottom squark 
mass which enter the corresponding e factors in the def- 
initions of Xb and Xk- As discussed in Sec. II A, such 
a splitting is compatible with the MFV ansatz for the 
squark spectrum as long as tan /3 is large enough that 
yb effects cannot be neglected. However, we find that 
even for a large mass splitting Xb — Xk holds, except 
for regions of parameter space with large and negative [i, 
such that e ( s )tan/3 ~ 0(— 1). Such regions of param- 
eter space are strongly constrained by perturbativity of 
the bottom Yukawa and vacuum stability considerations, 
as discussed in Sec. IV. If Xb — Xk, then the B — » tv 
decay gives stronger constraints than K — > fiv." 1 



4 For the special range 0.25 < Xg,Xx ;$ 0.30 and eotan/3 > 



In the following, we therefore concentrate on the con- 
straint from B — y tv on the MSSM parameter space. 
Apart from corners of parameter space with very large 
and negative eotan/3 < — 1, the charged Higgs contribu- 
tion interferes destructively with the SM (X B > 0), and 
leads to constraints in the My ± -tan /3 plane. These con- 
straints depend on other SUSY parameters only through 
the loop-induced tan /3 resummation factors e^ and are 
therefore robust in large parts of parameter space. 

The yellow lines in the left plot of Fig. 4 show the 
B — > tv constraints in the M^-tan /3 plane correspond- 
ing to the 3 choices of MSSM parameters (a), (b), and 
(c) given in Tab. I and already discussed in Sec. IIIB. 
For comparison, the constraint from direct searches in 
the t + t~ channel is also shown in gray. There are also a 
narrow strips of small Higgs masses and large tan /3 where 
the NP contribution to the B — > tv amplitude is twice as 
large as the SM contribution. This in turn implies that 
this region of parameters is in principle allowed by the ex- 
perimental data on5-> tv. It is in strong tension, how- 
ever, with the results from B — > Dtv, B — > D*tv, and 
K — > /iv and furthermore is excluded by direct searches. 

The dependence of the B —> tv constraints on the 
tan /3 resummation factors is stronger than the one of the 
direct searches in the t + t~ channel, especially for large 
values of tan j3. For large values of tan j3 and a positive 
(negative) value of eo the constraint can be weakened 
(strengthened) considerably. As e does not depend on 
A t , the constraint from B — > tv is to a large extent in- 
sensitive to the exact value of the light Higgs mass. Con- 
straints from direct MSSM Higgs searches are generically 
stronger for Ma < 800 GeV. While the latest results from 
direct MSSM Higgs searches in the t + t~ channel at the 
LHC end at Ma — 800 GeV, obviously no such restric- 
tion exists for the B — » tv constraints. Only very large 
values of tan/3 > 60, however, are typically probed by 
B — > tv for such large heavy Higgs masses. 

In corners of parameter space with very large negative 
fx, we can have eotan/3 < — 1 for values of tan/3 that 
are not extremely large and when the bottom Yukawa 
is perturbative [140]. Such a situation is shown in the 
right plot of Fig. 4, where [i = — 8 TeV and the resulting 
eo — —0.03. For tan/3 < 30, the charged Higgs still 
interferes destructively with the SM. For larger values of 
tan (3 > 30, however, the sign of X B flips, the interference 
becomes constructive, and the branching ratio is always 
enhanced. This behavior can be seen from the values of 
Rbtv indicated with the dotted contours in the right plot 
of Fig. 4. 

Note that vacuum stability requirements, however, 
strongly constrain very large and negative values of fi. As 
discussed in Sec. IV, we do not find viable regions of pa- 
rameter space where the bottom Yukawa has a negative 



— 1, the K — > fj,u constraint is stronger than B — > tv, but this 
parameter region is excluded by direct searches as discussed in 
the main text. 
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FIG. 4. Constraints in the Ma -tan f3 plane from the tree level B — > ri/ decay. The constraint from direct heavy Higgs searches 
is also shown in gray. The yellow solid, dotted and dashed contours in the left plot correspond to scenarios (a), (b), and (c) 
defined in Tab. I. The right plot shows a scenario with /i = —8 TeV, leading to a large negative eo such that the charged Higgs 
contribution interferes constructively with the SM in the region with tan/? > 30. The labeled contours indicate values for Rb T v 
Above the red horizontal line, the electroweak vacuum has a lifetime shorter than the age of the universe. 



sign with respect to the SM one, i.e. with e;,tan/3 < — 1. 
For B — > tv, the relevant parameter combination is 
eo tan /3. The horizontal red line in the right plot of Fig. 4 
marks the upper bound on tan (3 in the scenario with 
fj, = — 8 TeV, such that the electroweak vacuum remains 
stable on timescales of the age of the universe. Therefore, 
we see that eo tan (3 < — 1 is also excluded by vacuum sta- 
bility considerations. This conclusion holds beyond the 
discussed fi = — 8 TeV example. 



B. 



MA 4 



The B s — > fJ. n decay is a flavor changing neutral 
current process and correspondingly only induced at the 
loop level, both in the SM and the MSSM. In the SM, 
B s — > fi + fi~ is also helicity suppressed by the muon mass, 
resulting in a tiny SM prediction, at the level of 10~ 9 . 
Using the recently given precise value for the B s meson 
decay constant fg s = (227 ± 4) MeV [156] which is an 
average of several lattice determinations [153-155], and 
taking into account the effect of the large width difference 
in the B s meson system [173, 174], we have the branching 
ratio extracted from an untagged rate as [175] (see also 
[176]) 



BR{B S -> a* + aOsm = (3-32 ± 0.17) x 10" 



(41) 



Experimental searches for that decay have been carried 
out at DO [177] and CDF [178], and are ongoing at AT- 
LAS [179], CMS [180], andLHCb [80, 181]. Very recently, 



the LHCb collaboration reported first evidence for the 
B s — > /it + fi~ decay [80]. LHCb finds for the branching 
ratio the following value 



-9 



BR(B S -> m + aOc XP = (3.2 t\i ±S:I) x 10 



and gives the following two sided 95% C.L. bound 



1.1 x 10~ 9 < BR{B S -s- /z~V 



< 6.4 x 10" 



(42) 



(43) 



We use this bound in our analysis. Note that the upper 
bound in (43) is considerably weaker than the official 
combination of the previous LHCb result [181] with the 
ATLAS and CMS bounds [182]. 

For large values of tan/3, order of magnitude en- 
hancements of the BR(£> S — > /i + /i~) are possible in the 
MSSM [106, 183]. In the large tan/3 limit, the CP aver- 
aged branching ratio in the MFV MSSM can be written 
to a good approximation as 



Rb, 



BR(B s ^ f i+^-) 
BR{B S 



/""•"A* )SM 



|^| 2 + |1-^| 2 . (44) 



The MSSM contribution A is dominated by so-called 
Higgs penguins, i.e. the exchange of the heavy scalar 
H and pseudoscalar A with their 1-loop induced flavor 
changing b — > s couplings, that are parametrized by cfc 
given in (13). We find 



A = 



47r m| 



47r m- Bs e FC tp 



1 



(45) 



12 



The SM loop function Yq depends on the top mass and 
is approximately Yo ~ 0.96. Note that the MSSM con- 
tributions to B s — > do not decouple with the scale 
of the SUSY particles, but with the masses of the heavy 
scalar and pseudoscalar Higgs bosons Mfj ~ M\. Due to 
the strong enhancement by tan 3 j3, the large tan (3 regime 
of the MSSM is highly constrained by the current exper- 
imental results on BR(£> S — > We remark, how- 
ever, that cfc in the numerator of (45) is a sum of several 
terms (see (23)) each of which depend strongly on several 
MSSM parameters. In addition, cancellations among the 
different terms can occur in certain regions of parameter 
space, rendering the B s — > /i+ 'fi~ constraint very model 
dependent, even in the restrictive framework of MFV. 
Additional contributions to B s — > can arise from 
charged Higgs loops [187]. They interfere destructively 
with the SM contribution and scale as (tan/3) 2 /Af^ ± . 
Typically, their effect is considerably smaller compared 
to the SUSY contribution in (45). 

We stress that there is a simple mathematical lower 
bound of Rb s hh = 1/2 in (44) that is saturated for 
A — 1/2. In this case, the SUSY contribution partially 
cancels the SM amplitude, but simultaneously generates 
a non-interfering piece that cannot be canceled. This 
lower limit provides a significant threshold for experi- 
ments searching for BR(£> S — > /i + /i~): not only is the 
SM branching fraction a meaningful value to test experi- 
mentally, but the potential observation of the branching 
fraction below one half of the SM value would strongly 
indicate NP and imply departure from the MSSM with 
MFV. Note that the current 2a lower bound from LHCb 
on the branching ratio is below 1/2 of the SM value and 
therefore does not lead to constraints in our framework, 
yet. 

In Fig. 5, we show the constraints from B s — > in 
the M^-tan/3 plane. The red solid, dotted and dashed 
contours correspond to scenarios (a), (b), and (c) of 
Tab. I. The dash-dotted contour corresponds to scenario 
(d), with all MSSM parameters as for the solid con- 
tour, but with a negative sign for the trilinear coupling. 
For comparison, the constraints from direct searches are 
again shown in gray. As expected, we observe a very 
strong dependence of the B s — > /J, + fJ,~ bounds on the 
choices of the remaining MSSM parameters, particularly 
the sign of \iA t . Note that in the considered scenarios, 
we assume degenerate squarks such that the only term 
entering epc is from the irreducible Higgsino loop contri- 
bution, , whose sign is dictated by iiA t . For positive 
(negative) \j,A t the NP contribution interferes destruc- 
tively (constructively) with the SM amplitude. Since the 
lower bound on BR(£? S — » n + ^jT) from LHCb is still be- 
low half of the SM value, destructively interfering NP is 
much less constrained than constructively interfering NP. 

The plots of Fig. 6 show in red the constraints from 
B s — > in the plane of the third generation squark 

masses and the Higgsino mass parameter fj,. The gray 
horizontal band corresponds to the constraint from di- 
rect searches of charginos at LEP that exclude \fi\ < 
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FIG. 5. Constraints in the Ma -tan /3 plane from the B s — > 
fi + fi~ decay. The red solid, dotted, dashed and dash-dotted 
contours correspond to scenarios (a), (b), (c) and (d), as de- 
scribed in the text. The gray region is excluded by direct 
searches of MSSM Higgs bosons in the H/A — > t + t~ chan- 
nel. 



100 GeV [184, 185]. In these plots, we fix M A = 800 GeV, 
tan j3 = 45 (fully compatible with the B — > tv constraint 
and not yet constrained by direct searches) , and gaugino 
masses with 6Mi = 3M 2 = M 3 = 1.5 TeV. As in all the 
other plots, we vary the trilinear couplings A t = Ab = A T 
throughout the plot such that the lightest Higgs mass is 
Mh = 125 GeV. The values for A t are indicated in the 
plots by the vertical dotted contours. The two plots cor- 
respond to positive and negative values of the ^4-terms. 
In the gray region in the lower left corners of the plots, the 
sbottom loop corrections to the lightest Higgs mass be- 
come so large that the lightest Higgs mass is always below 
Mh < 125 GeV for any value of A t , taking into account 
a 3 GeV theory uncertainty. We checked that varying 
the light Higgs mass between 122 GeV < M h < 128 GeV 
can change the values of A t by around 25% in each di- 
rection and therefore can affect the constraints derived 
from B s — > fi + pT at a quantitative level. However, the 
qualitative picture of the constraints and the interplay 
of the SUSY contributions to B s — > p + n~ , as discussed 
below, are unaffected by this variation. 

The solid contours are obtained under the assumption 
that the masses of the first two generation squarks are 
equal to the third generation, while for the dashed and 
dotted contours we assume the first two generations to 
be heavier by 50%. For the dashed contours, we as- 
sume the splitting for the left-handed squarks to be fully 
aligned in the up-sector, such that gaugino-squark loops 
also contribute to epc with £ = 1 (see (23) and (25)). 
We set C = 0.5 for the dotted contours, such that only 
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FIG. 6. Constraints in the mq 3 -(i plane from the B s -> [i + fi~ decay, with fixed M 3 = 3M 2 = 6M1 = 1.5 TeV, M A = 800 GeV 
and tan/3 = 45. The solid bounded regions correspond to a degenerate squark spectrum. The dashed and dotted bounded 
regions correspond to choosing the first two squark generations 50% heavier than the third generation squark masses, with an 
alignment of £ = 1 and £ = 0.5, respectively. The gray horizontal band corresponds to the constraint from direct searches of 
charginos at LEP. The vertical dotted lines show contours of constant At such that Mh = 125 GeV. In the gray regions in the 
lower left corners, the lightest Higgs mass is always below Mh < 125 GeV, taking into account a 3 GeV theory uncertainty. 



half of the squark mass splitting induces flavor viola- 
tion in the down-sector. For negative A t , the obtained 
bounds show a strong dependence on the value of C- The 
BR(£> S — > /i + /i~) bounds in Fig. 6 clearly display the 
non-decoupling behavior mentioned above. Due to this 
non-decoupling, the BR(.B,; — > results can con- 

strain SUSY parameter space in regions that are beyond 
the current and expected future reach of direct searches. 

A crucial element of our analysis is the viability of the 
cancellation of the SUSY contribution to the B s — > 
branching ratio. This cancellation is driven by the pres- 
ence of 6fc m (45), which is schematically given in (23) 
and its various contributions are detailed in (21), (25) 
and (26). First, in the following discussion, we neglect 
the wino contribution given by (26), which is generally 
smaller than the gluino contribution. This is due to the 
smallness of M 2 and a in (26) compared to M 3 and a s 
in (25) (of course, our numerical analysis always includes 
the wino contribution). Since each SUSY contribution 
is proportional to (i, we see that switching the sign of 
fi changes the relative sign between the SUSY and SM 
amplitudes. Furthermore, by switching the sign of A tl 
between the left and right panels of Fig. (6), we change 
the relative sign between the gluino contribution and the 
Higgsino contribution. Thus, for a particular choice of 
sign(A t ) and sign(/x), we can exploit a cancellation be- 
tween the gluino vs. Higgsino loop, diminishing the mag- 
nitude of the SUSY contribution, and a second cancella- 
tion between the overall SUSY contribution and the SM 
amplitude. In particular, even if the magnitude of the 



SUSY contribution is by itself larger than the SM con- 
tribution, we can exercise the second cancellation where 
the SUSY amplitude overshoots the SM one. 

These cancellations are clearly in effect in the left and 
right panels of Fig. 6. We first focus on the regions 
bounded by solid lines, which correspond to degenerate 
squark masses. This implies that the SUSY contribution 
dominantly arises from in (21). In the upper half of 
the left panel corresponding to positive A t and positive /x, 
the SUSY contribution cancels with the SM contribution 
and always leads to a BR(B S — > /i + /i~) below the current 
bound. In the lower half of the left panel, with positive 
A t and negative [i, the Higgsino contribution adds con- 
structively with the SM contribution, leading to signifi- 
cant constraints. In the upper half of the right panel, the 
Higgsino contribution also adds constructively with the 
SM, leading again to a bound. This bound is less strin- 
gent compared to the positive A t and negative [i case, 
because for positive fi, the e;, and 6q terms in (45) lead 
to a suppression of the SUSY amplitude. Finally, in the 
lower half of the right panel, with negative A t and nega- 
tive ii, the Higgsino contribution interferes destructively 
with the SM. The constraint is non-vanishing, however, 
because for negative /j, 7 the tan /3 resummation factors, 
given in (45), enhance the SUSY amplitude such that it 
can be more than twice as large as the SM amplitude. 

When we include squark splitting, we further 
strengthen the SUSY contribution for positive A t , be- 
cause the gluino and Higgsino contributions add con- 
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FIG. 7. The Myt-tan ft plane in view of projected constraints 
from the BR(_B S — > fj, + fi~), assuming a future ±0.5 x 10~ 9 
uncertainty in the measurement with the SM prediction as 
the central value. The green shaded regions between and 
below the solid and dashed contours correspond to values for 
tan/? and Ma allowed in scenarios (a) and (e), as defined in 
Tab. I. The gray region is excluded by current direct searches 
of MSSM Higgs bosons in the H/A -+ t+t~ channel. 



structively. Hence the overall SUSY+SM interference is 
more restricted. The bounds due to this splitting in the 
masses are shown by the regions enclosed by the dashed 
and dotted lines in Fig. (6). For negative A t , shown in 
the right panel, the gluino contribution partially cancels 
the Higgsino contribution, leading to a weaker constraint. 
The effect of the gluino contributions decreases for larger 
gluino mass, M 3 . 

In tandem, the complementary views provided by the 
different panels of Figs. 5 and 6 clearly demonstrate that 
certain choices of SUSY parameters relax the constraints 
considerably. For example, with Ma — 800 GeV and 
tan/3 = 45, the region of parameter space with positive /j, 
and positive A t is robustly unconstrained from the B s — > 
limit. Moving from top to bottom along a constant 
A t contour in the left plot of Fig. 6 corresponds to a rapid 
coverage of the tan /3 vs. Ma plane from the (b) to (a) 
to (c) exclusion regions. 

Regions of parameter space with destructive interfer- 
ence between SM and SUSY amplitudes (i.e. the regions 
with positive fxAt) will be constrained significantly if a 
lower bound of BR(B S — > /i + /i~) above half of the SM 
prediction is established in the future. We illustrate this 
in the plots of Figs. 7 and 8, which assume a measurement 
of BR(B S — > n + fi~) at the SM expectation as a central 
value with an experimental uncertainty of ±0.5 x 10~ 9 . 
Such a precision is expected to be achieved by LHCb at 
the end of the 13 TeV run with a combined analysis of 



1 ftr 1 of 7 TeV data, 1.5 ftr 1 of 8 TeV data, and 4 
fb" 1 of 13 TeV data [186]. The plots in Figs. 7 and 8 
show in green the regions in the M^j-tan/3 and TnQ 3 —ji 
planes that are allowed by the expected results on the 
B s — > decay As shown in Fig. 7, apart from the 

allowed regions with large Ma and small tan j3, there are 
also strips with large Ma and large tan /3 where the ex- 
pected bounds from B s — > fJ- + fJ-~ can be avoided. In 
these regions, the SUSY amplitude has approximately 
the same size as the SM amplitude but is opposite in 
sign. According to (44), this leads to a branching ratio 
close to the SM prediction. 

For the example parameter point with Ma — 800 
GeV and tan/3 = 45, the projected lower bound on 
BR(B S — » leads to very strong constraints in the 

m Q3~ A* plane for positive /iA t . Indeed, for Ma = 800 
GeV and tan/3 = 45, and given the assumed exper- 
imental precision, charged Higgs loop contributions to 
B s — > already lead to a non-negligible suppres- 

sion [187], leaving hardly any room for destructively in- 
terfering SUSY contributions. Only if the SUSY contri- 
bution is so large that A ~ — 1 does the parameter space 
open up again. The corresponding regions that are ex- 
cluded by the assumed lower bound are clearly visible in 
the white region of the upper half of the left plot and the 
upper white region in the lower half of the right plot in 
Fig. 8. 



C. B — >■ X s ~f 

The loop induced B — > X s j decay is also highly sen- 
sitive to NP effects coming from SUSY particles. The 
NNLO SM prediction for the branching ratio reads [188] 
(see also [189, 190]) 



BR(B -> Y s7 ) S m = (3.15 ± 0.23) x 10" 



(46) 



On the experimental side, BaBar recently presented up- 
dated results for the branching ratio [83]. Including this, 
the new world average reads [191] 



BR(B -+ Y s7 ) OX p = (3.43 ± 0.22) x 10" 



(47) 



which is slightly lower than the previous world average 
and is in very good agreement with the SM prediction. 
In the MSSM with minimal flavor violation and no new 
sources of CP violation, the branching ratio can be writ- 
ten as [192] 



Rbsj — 



BR{B -> X sl ) 
BR{B -+ X sl ) 



SM 



NP 



0.61 C 8 NP + 0.74C 7 NP C 8 NP 



~ 1 - 2.55 Cy 

±1.57 (C 7 NP ) 2 ± 0.11 (C* 8 NP ) 2 , (48) 

where C^g are the NP contributions to the magnetic and 
chromo-magnetic b — > sj operators evaluated at the scale 
160 GeV. 
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FIG. 8. The mQ 3 -\i plane in view of projected constraints from the BR(_B S — > fi + fj,~), assuming a future ±0.5 x 10~ 9 uncertainty 
in the measurement with the SM prediction as the central value. We fixed M3 = 3M2 = 6M1 — 1.5 TeV, Ma = 800 GeV 
and tan /3 = 45. The green shaded regions between the solid contours correspond to values for rriQ 3 and fi allowed for a 
degenerate squark spectrum. The green shaded regions between and above the dashed contours are allowed if the first two 
squark generations are 50% heavier than the third generation squark masses, with an alignment of £ = 1. The gray horizontal 
band corresponds to the constraint from direct searches of charginos at LEP. The vertical dotted lines show contours of constant 
At such that Mh = 125 GeV. In the gray regions in the lower left corners, the lightest Higgs mass is always below Mu < 125 GeV, 
taking into account a 3 GeV theory uncertainty. 



Apart from the B — > X s j decay, the modifications of 
the Wilson coefficients C7 and Cg, also enter predictions 
of observables in the B -+ !<*£+£- decay. In our MSSM 
setup with minimal flavor and CP violation, we only have 
real NP contributions to C7 and C%. In this framework, 
the experimental data on B -+ K*£ + £~ does not put 
additional restrictions, once the bounds from BR(_B — > 
X s -f) are taken into account [175, 193]. Therefore, we 
focus only on the B — » X s ^f decay. 



The SUSY contributions to Cf g come from charged 
Higgs-top loops, neutral Higgs-bottom loops, Higgsino- 
stop loops, and gaugino-squark loops. As with the 
Higgs-fermion couplings, we take into account the most 
generic MFV structure of the squark masses and consis- 
tently consider splittings between the first two and the 
third generation squarks in the left-handed as well as the 
right handed sector. The resulting dominant MSSM con- 
tributions to Cr,s read 



°7,8 



°7,8 



r<9 

i/3 



°7,8 



1-^ 

l + e b tp (l + e b tp){l + e t p ) I 2M 
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M^M 3 fl 8 (m 2 Q3 ,m 2 D3 ,M^) 
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(49) 
(50) 

(51) 
(52) 
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FIG. 9. Constraints in the Ma -tan /3 plane from the B — > X s ^y 
decay. The orange solid, dotted, dashed, and dash-dotted 
contours correspond to scenarios (a), (b), (c), and (d) as de- 
scribed in the text. The gray region is excluded by direct 
searches of MSSM Higgs bosons in the H/A — > t + t~ chan- 
nel. 



The first term in (49) corresponds to contributions from 
a charged Higgs loop. The loop functions, hr,s depend 
on the ratio of the top mass and the charged Higgs 
mass, r t — ml /Mtt± , and for r t — 1 are given by 
h 7 (l) = -7/18 and h 8 (l) = -1/3. Their full analytical 
expressions can be found in the appendix. The second 
term in (49) arises from neutral heavy Higgs loops. It 
is strongly suppressed by the bottom quark mass and is 
only important for very large tan/3. The loop functions, 
depend on the ratio of the bottom mass and the 
charged Higgs mass and since m^/M^ ± <C 1, they are 
very well approximated by zj — — jg and zg = g . 



Contributions from Higgsino-stop, gluino-down 
squark, and Wino-down squark loops are shown 
in (50), (51), and (52), respectively. We do not write the 
typically negligible bino contributions. 

For a degenerate SUSY spectrum with mass m, the 
loop functions entering the Higgsino and gaugino contri- 
butions reduce to 

2 1 

27?^ 
5 1 
18^4 



ft 



5 1 

36 m 4 
1 1 

12^ : 



/I 



h 



7 1 

24^4 
1 1 

"8^4 ' 



Their full analytical expressions are collected in the ap- 
pendix. In contrast to the Higgs penguin contributions 
to B s — > fi + fi—, the SUSY loop contributions to b — > S7 
do decouple with the SUSY scale. 

The first terms in (51) and (52) correspond to 1-loop 
flavor changing gaugino contributions. They vanish for 



m Q'i = m Qi *- e - if there is no splitting between the 
first two and the third generations of left-handed squark 
masses. In the presence of a splitting, the parameter 
C again parametrizes the alignment of the left-handed 
squark mass matrix. As mentioned before, if the splitting 
is generated by RGE running we expect 1/2 < ( < 1. 
The second terms in (51) and (52) are formally 2-loop 
contributions but they can be relevant for large tan /3. 
They do not vanish for degenerate masses [110, 143]. 

Similarly to B s — > /i + /J~, the MSSM contribution to 
B — > Y s 7 is a sum of several terms that depend sensi- 
tively on many parameters, particularly the signs of fi 
and A t . 

In Fig. 9, we show in orange the constraints from 
B — > X s j in the M^-tan /3 plane obtained analogous 
to the B s — > fi + (i~ constraints discussed previously. The 
plots of Fig. 10 show the B — > X s -f constraints in the 
plane of the third generation squark masses and the Hig- 
gsino mass parameter fj,, again in complete analogy to 
the B s — > constraints. 

We can again see the connection between the con- 
straints in the tan/3 vs. Ma plane, given in Fig. 9, and 
the 11 vs. rriQ 3 plane, given in Fig. 10. The squark masses 
are fixed to 2 TeV in Fig. 9. This causes the stop-chargino 
contribution to be essentially negligible, and hence we 
are only constrained by the Higgs contribution in the 
low Ma and large tan/3 regions. For heavy squarks and 
low tan/3, the bound on the charged Higgs mass is ap- 
proximately independent of the other SUSY parameters 
and is given by M H ± > 300 GeV. For large tan /3, the 
rcsummation factors in (49) become relevant. The most 
important effect arises from the factors e' and £{, in the 
first term in (49). For negative /z, e' Q and eb are nega- 
tive and therefore the bounds become stronger for larger 
tan /3 in scenario (c) . For positive [i (scenarios a, b, and 
d) instead, the bounds are relaxed for large tan/3. As the 
dominant gluino contribution to e' and tb grows with fi 
the B — > X s j constraint is weakest in scenario (b) that 
has the largest [i — 4 TeV. For the heavy squark masses 
chosen in Fig. 9, the direct searches for MSSM Higgs 
bosons give stronger constraints compared to B — > X s j 
except for small values of tan /3. 

In the plots of Fig. 10, the variation of the squark 
masses allows the stop-chargino contribution to become 
important for small rriQ 3 , demonstrating that the tan /3 
vs. Ma projection insufficiently illustrates the B — > X s ^ 
constraint. Partial cancellations are again in effect, and 
we describe the relative signs of the various contributions 
in the following. Apart from extreme regions of parame- 
ter space, the charged Higgs contribution interferes con- 
structively with the SM and enhances BR(£? — > X s j). 
However, for the case shown, Ma = 800 GeV, this con- 
tribution is small. For positive (negative) (nA t ), the Hig- 
gsino loop contribution come with same (opposite) sign 
with respect to the SM. Among the gaugino contribu- 
tions, the dominant one is typically the 1-loop gluino con- 
tribution. If a splitting in the left-handed squark masses 
is induced radiatively, its sign depends, for positive M3, 
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FIG. 10. Constraints in the rriQ 3 -(j, plane from the B — > X s ~{ decay, for fixed M3 = 3M2 = 6M1 = 1.5 TeV. The solid bounded 
regions correspond to a degenerate squark spectrum. The dashed and dotted bounded regions correspond to choosing the 
first two squark generations 50% heavier than the third generation squark masses, with an alignment of f = 1 and £ = 0.5, 
respectively. The gray horizontal band corresponds to the constraint from direct searches of charginos at LEP. The vertical 
dotted lines show contours of constant At such that Mh = 125 GeV. In the gray regions in the lower left corners the lightest 
Higgs mass is always below Mh < 125 GeV, taking into account a 3 GeV theory uncertainty. 



only on the sign of /i. For positive (negative) /x, gluinos 
interfere destructively (constructively) with the SM. 

The plots of Fig. 10 clearly show the decoupling be- 
havior of the MSSM contributions to the b — > 57 transi- 
tion. For a degenerate squark spectrum (m.Q 3 — ttiq = 
mu a = nijj = niflj = rriD = rh) and a heavy charged 
Higgs, the bound from BK(B — > X s ^) hardly constrains 
the MSSM parameter space beyond squark masses that 
are already excluded by direct SUSY searches, namely 
™ ^ C(l TeV). In the presence of a mass splitting be- 
tween the first two and the third generations of squarks, 
the B — > X s "/ constraint can become relevant for nega- 
tive A t , since the gluino and Higgsino contributions add 
constructively. Squark masses significantly above 1 TeV 
can be probed in that case. For positive values of A t , on 
the other hand, the gluino and Higgsino loops partially 
cancel and the bound from B —¥ X s j is barely relevant. 



D. Discussion of RGE Effects 

Our phenomenological analysis of MSSM mass param- 
eters serves our purpose of understanding the flavor con- 
straints on the low energy MSSM spectrum. However, 
we also want to connect these constraints to parame- 
ters of a high scale SUSY parameter space. To this 
end, we consider a typical example in the large tan/3 
and Ma region compatible with direct H/A — > r + r~ 
searches at the LHC. We show typical mass differences 



between soft parameters for squarks in the plane of the 
mSUGRA boundary conditions, mo and m^, fixing the 
remaining mSUGRA parameters to Aq — ±2 TeV and 
tan/3 = 45. We also chose the SUSY breaking scale to be 
the GUT scale, 10 16 GeV. We deviate slightly from the 
strict mSUGRA prescription and work in a non-universal 
Higgs mass (NUHM) scenario by fixing the Higgs soft 
mass m 2 H = 1 TeV 2 and adjusting rrijj at the high 
scale to obtain Ma within 10% of 800 GeV at the low 
scale of Q = 1 TeV. Using these boundary conditions and 
the usual low energy Yukawa constraints derived from 
fermion masses run to Q = 1 TeV, we numerically solve 
the RGE system dictated by 2-loop running from [194] 
and 1-loop radiative corrections from [195]. Our choice 
of Ao typically gives the lightest SM-like Higgs a mass of 
122± 2 GeV. For the bulk of the region in the (mo, m 1( / 2 ) 
plane, adjusting Aq (in particular, A t ) to obtain a Higgs 
mass of 125 GeV changes the quantitative picture by less 
than a few percent. For very small mo and mi/2, how- 
ever, where some squarks or sleptons become close to 
tachyonic, the mass splittings can vary significantly as 
result of changing A . 

We highlight that the B observable constraints can 
vary significantly as a result of Yukawa-induced squark 
mass splittings inherent in RG running, as seen in Fig. 6 
and Fig. 10, respectively. In particular, the most signifi- 
cant mass splittings among the squarks occur as a result 
of the top and bottom Yukawas, where a significant en- 
hancement of the bottom Yukawa occurs for large tan /3. 

We can obtain a semi-analytic understanding of the re- 
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FIG. 11. Contours of AQ13 (top left), AC/13 (top middle), AD13 (top right), AQC/3 (bottom left), and AQD 3 (bottom right) 
in the {m ,m 1/2 ) plane, fixing A = 2 TeV, (m° H J 2 = 1 TeV 2 , tan/3 = 45, and requiring M A = 800 GeV within 10%. 



suiting mass splittings, following the simplified 1-loop RG 
analysis of [196-198]. Neglecting the first- and second- 
generation Yukawa couplings and a 2 contributions, we 
have 

m Q» W ~ m Qs (°) + 7 «3 + I a2 -It -h, (53) 
<(t)»<(0)+^-W tl (54) 



(55) 



and 
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where t = corresponds to the GUT scale. The anal- 



ogous m Ql , 



and mjj i approximations can be ob- 



tained from the above by neglecting the I t and If, contri- 
butions. 

For the trilinear couplings, neglecting a\ and A T , 
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The most relevant mass splittings for our analysis are 
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From these relations we see that AC/13 + AD13 ~ 2AQi 3 , 
where AZ?i 3 is small compared to AC7i 3 for small tan /3. 
We also expect AQC7 3 = —AQD 3 for small m^- These 
relations for the various mass splittings, based on 1-loop 
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FIG. 12. Same as Fig. 11 except A = -2 TeV. 



semi- analytic results, are borne out in our numerical re- 
sults, which are calculated from 2-loop RG running, and 
are shown in Fig. 11 and Fig. 12. 

For a gluino with mass M 3 = 1.5 TeV as we consid- 
ered in the previous sections, we have m 1 / 2 — 500 GeV 
and therefore a significant splitting is induced among the 
squark masses in running down from the GUT scale. The 
most important splitting in the general discussion of the 
previous sections is AQ13, as it leads to gaugino loop 
contributions to FCNCs. In our mSUGRA setup, it is 
typically around 25% for positive A and 35% for nega- 
tive A . The splitting between the masses of the squarks 
decreases for larger mi/ 2 . This is due to the universal 
SU(3) contribution, I a3 , to m^, ttzq, mjj 3 , rnu, n%D 3 , 
and rri£), which dominates for large m-i/2- 

From the approximate expressions above, we can also 
estimate the size of £ resulting from running. We have 



c 



It + h 



(67) 



Even though we chose a large value of tan/3 = 45 for 
the examples shown, the bottom Yukawa effects are lim- 
ited. Note that for the parameter region explored here, 
£ ~ 80%, which means that the squark mass splitting 
is dominantly driven by the top Yukawa and therefore 
aligned in the up-sector. For smaller tan /3, the align- 
ment parameter £ is even closer to 1. 

Note that the gaugino loop contributions to FCNCs 



depend approximately on the product £ x AQ13. In the 
mSUGRA scenario discussed here, we find to a good ap- 
proximation £ x AQ13 ~ AC/13/2. In more generic setups 
however, this relation does not hold and we will continue 
to discuss the gaugino loop contributions to FCNCs in 
terms of £ and AQ13 separately. 

In the plots of Fig. 13, we show again the constraints 
from B s — > fi + pr and B — > X s ^ in the mg 3 -/i plane, 
this time setting the various mass splittings according to 
our results of the mSUGRA RGE running. In particu- 
lar, we use AQ13 = 0.35, AC/13 = 0.6, AL>i 3 = 0.15, 
AQC/3 = 0.35, and AQD3 = —0.25, which are typical 
values for — 500 GeV and negative A t . As we saw 
in the previous sections for positive At, the B s — > 
constraint depends very mildly on the squark mass split- 
ting and the B — > X s j constraint is barely relevant. 
Therefore, we restrict ourselves to negative A t . For com- 
parison, the solid contours indicate again the constraints 
obtained for a degenerate squark spectrum. The dot- 
ted contours corresponds to keeping all third generation 
squarks degenerate and only implementing the splitting 
between the first two and the third generation as given 
by the RGE running. The dashed contours correspond 
to the situation where all squark mass splittings are as 
dictated by the RGE running. The former case behaves 
as expected given the analysis of Sees. VB and VC. For 
the latter case, however, once mass splittings between 
the different types of third generation squarks are also 
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negative A t negative A, 




FIG. 13. Constraints in the rriQ 3 -fi plane from B s — > M + M (left) and B — > X s ^ (right). The solid bounded regions correspond 
to a degenerate squark spectrum. The dashed and dotted bounded regions correspond to mass splittings in the squark spectrum 
implied by RGE running. In particular AQ 13 = 0.35, A(7i 3 = 0.6, A£>i 3 = 0.15, AQU 3 = 0.35, and AQD 3 = -0.25 for the 
dashed contours and AQ13 = 0.35, AU13 = 0.6, AD13 = 0.15, and AQU3 = AQD3 = for the dotted contours. The gray 
horizontal band corresponds to the constraint from direct chargino searches. The vertical dotted lines show contours of constant 
At such that Mh = 125 GeV. In the gray regions in the lower left corners the lightest Higgs mass is always below < 125 GeV, 
taking into account a 3 GeV theory uncertainty. 



considered, an additional effect arises. As can be seen 
from (54) and confirmed in the lower left plots of Figs. 11 
and 12, the right handed stop is typically significantly 
lighter than the other third generation squarks. The light 
right-handed stop then increases the chargino-stop loop 
contributions to B s — > n + fJ-~ and to B — > X s j leading 
overall to stronger constraints compared to the case of 
degenerate third generation squarks. 



Two of the most important quantities dictated by 
RGEs for flavor observables are the values of AQ 13 and 
£. Within the assumption of flavor universality at the 
messanger scale, AQ13 and £ depend mainly on the mes- 
senger scale, tan /? and the ratio of gluino mass to squark 
masses. Lowering the messenger scale from the GUT 
scale as well as increasing the gluino mass decreases the 
splitting AQ13, but leaves C approximately invariant. 
Smaller (larger) values of tan /? would decrease (increase) 
AQ13 and simultaneously increase (decrease) £, leaving 
the product £ x AQ13 — AU13/2 approximately invari- 
ant. As we saw, making the splitting smaller strength- 
ens the BR(B S — > /i + /i~) constraint for negative A t , but 
increasing £ will relax it. The effect of these two quan- 
tities is exactly opposite on the constraints coming from 
BR(B — > X s j). This complimentary behavior implies 
that even varying the messenger scale and tan j3, these 
two flavor observables will be able to constrain the pa- 
rameter space efficiently. 



VI. DARK MATTER DIRECT DETECTION 

The lightest neutralino in the MSSM is an excellent 
thermal dark matter candidate. The lightest neutralino is 
a weakly interacting massive particle (WIMP) and there- 
fore generically leads to roughly the correct order of mag- 
nitude for the observed dark matter relic density. This 
is particularly true in the well-tempered neutralino sce- 
nario [199], where the lightest neutralino is a mixture 
of the Bino and Wino or the Higgsino. In the follow- 
ing we do not assume any specific mechanism by which 
the correct dark matter relic abundance is achieved, but 
simply assume that the lightest neutralino in the MSSM 
accounts for the dark matter in the universe [200]. 

Neutralinos interact with SM matter and therefore 
dark matter direct detection limits can be used to put 
bounds on the MSSM parameter space, complementary 
to the bounds from direct searches and low energy flavor 
observables [59, 68, 201-209]. 

The XenonlOO Collaboration recently set very strin- 
gent limits on the spin-independent elastic dark matter 
nucleon scattering cross section [84, 211]. For dark mat- 
ter masses of O(100 GeV), the bounds are as strong as 
<7 < 2 x 10 _45 cm 2 , assuming canonical values for the lo- 
cal dark matter density, the local circular velocity and 
the Galactic escape velocity. Interpreted in the context 
of the MSSM with neutralino dark matter, these bounds 
are starting to probe significant parts of the parameter 
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FIG. 14. Constraints in the Ma -tan /3 plane from Dark Matter direct detection. The green solid, dashed, dotted and dash- 
dotted contours correspond to different values of /x as indicated. The gray region is excluded by direct searches of MSSM Higgs 
bosons in the H/A — > t + t~ channel. 



space. 

The spin-independent elastic neutralino-proton cross- 
section can be written as 



7r(M x + TO p ) ; 



-f 



(68) 



where M x is the mass of the lightest neutralino, m p is 
the proton mass, and 



fT g c q + 97 /t 



\ " 



q=c,b,t 



(69) 



The non-perturbative parameters and f^ G = 1 — 

/t — /t — /t come from the evaluation of nuclear ma- 
trix elements. We use the latest lattice determinations 
in our numerical analysis [210] 



= /£ d = 0.028 , /£ s = 0.0689 



(70) 



These values are expected to be affected by considerable 
uncertainties. We assume isospin symmetry when apply- 
ing the XenonlOO bounds. 

For large tan/3, the dominant contributions to the co- 
efficients, c qi parametrizing the neutralino-quark cou- 
plings, typically come from the i-channel exchange of the 
heavy scalar H and read 



c d — c s 



9l 



ti- 



lt 



AM% 1 + e s t Ml - ' 
U4 2 H 1 + e b t p Ml - /i 2 • 



(71) 
(72) 



The i-channel exchange of the SM-like Higgs affects all 
c q approximately equally: 



Mi 



"« ~ AMI M l 



(73) 



While the are not enhanced by tan /3, bounds on 
the direct detection cross section have become so strong 
that the i-channel exchange of the SM-like Higgs is also 
probed. 

The above expressions hold in the large tan /3 limit and 
assume the lightest supersymmetric particle to be mainly 
a bino-higgsino mixture with M\ ^= fi. In our numerical 
analysis, we go beyond the large tan/3 limit: we work 
with neutralino mass eigenstates and include the effects 
from s-channel squark exchange, though these are always 
very suppressed by the squark masses. 

As is evident from (71), (72), and (73), the neutralino- 
proton cross section depends strongly on Mi and /i. This 
can be also seen from the plots of Fig. 14, which show 
in green the regions in the canonical M^-tan/3 plane 
that are excluded by the XenonlOO constraints. In the 
left plot, the bino mass is set to Mi = 100 GeV with 
M 2 = 2A/i and the solid, dashed, dotted and dash-dotted 
contours correspond to (i = 250 GeV, 300 GeV, 400 GeV, 
and 600 GeV, as indicated in the plot. In the right plot, 
we choose a larger bino mass of Mi = 600 GeV, with 
M 2 = 2Mi again, and the solid, dashed, dotted and dash- 
dotted contours correspond to ^ = 250 GeV, 500 GeV, 
750 GeV, and 1000 GeV. Dependence on other SUSY 
parameters enters at the loop level through the £j fac- 
tors in (71) and (72) and is therefore very moderate. In 
these plots we fix a common squarks mass, rh = 2 TeV, 
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Mi (TeV) 

FIG. 15. Constraints in the Mi—fi plane from Dark Matter 
direct detection. The solid and dashed contours correspond to 
different choices for Ma and tan /3 as denned in the text. The 
horizontal gray band is excluded by direct chargino searches. 



a gluino mass of M3 = 1.5 TeV and A t = A\, = A T such 
that the lightest Higgs mass 125 GeV. 

The strongest constraints arise if binos and higgsi- 
nos are maximally mixed, i.e. for Mi ~ fi. Indeed, if 
Mi = [i, we find, independent of the values of Ma and 
tan /3, that the exchange of the SM-like Higgs leads to di- 
rect detection cross sections that are already ruled out by 
the current bounds in the full range of neutralino masses 
up to 1 TeV. Away from bino-higgsino degeneracy, re- 
gions of parameter space open up. Still, for small heavy 
Higgs masses and large values of tan /3, the heavy Higgs 
exchange contributions can be sizable and lead to impor- 
tant constraints in the M^-tan /3 plane, as long as \x and 
Mi < 1 TeV. In the excluded regions with small tan/3 
and a large heavy Higgs mass, the constraint arises from 
the exchange of the light Higgs. 

The plot of Fig. 15 shows the direct detection con- 
straints in the Mi~\x plane for 2 different points in the 
Ma -tan /3 plane. The solid, dashed and dotted contours 
correspond to Ma = 800 GeV and tan/3 = 45 (scenario 
i) and Ma = 1 TeV and tan/3 = 10 (scenario ii), both 
compatible with current direct searches. As already men- 
tioned, the strongest constraints arise along the Mi ~ fi 
line. Interestingly, the constraints for negative values of 
/1 are considerably weaker, because for negative /1, the 
heavy Higgs and SM like Higgs contributions interfere 
destructively. Observe that this behavior is opposite to 
that of the constraints coming from B s — > and 
B — > AT S 7, which are currently weaker for positive fi (and 
positive A t ). 

Note that the bounds from dark matter direct detec- 
tion not only depend very strongly on various MSSM pa- 



rameters, but are also affected by various uncertainties, 
e.g. from the nuclear matrix elements, and astrophysical 
uncertainties, in particular the dark matter velocity dis- 
tribution. Moreover, they also depend crucially on the 
assumption that the dark matter of the universe indeed 
consist entirely of MSSM neutralinos. If neutralinos only 
make up a (small) fraction of the dark matter, the bounds 
can be relaxed considerably and even avoided completely. 



VII. CONCLUSIONS 

In this work, we evaluated the status of the minimal 
supersymmetric standard model with minimal flavor vi- 
olation in light of the recent Higgs discovery as well as 
constraints from collider searches, flavor measurements, 
and dark matter direct detection experiments. In con- 
cert, these complementary probes provide valuable con- 
straints on the MSSM parameter space. In particular, we 
showed that flavor bounds can be stronger than bounds 
from direct searches for heavy MSSM Higgs particles or 
supersymmetric particles, even in the restrictive frame- 
work of MFV. 

Throughout our analysis, we consistently implemented 
the most general structure of the soft SUSY breaking 
terms compatible with the MFV ansatz, i.e. allowing 
splitting between the first two and the third generations 
of squarks. We demonstrated that, in addition to the 
typical pMSSM parameters, an additional parameter, 
reflective of the alignment of the mass splitting of the 
left-handed squarks, is required to discuss the flavor phe- 
nomenology of this framework. In the presence of such 
splitting, this parameter controls the size of gaugino- 
squark loop contributions to FCNCs. Possible cancel- 
lations between gaugino and higgsino loop contributions 
have a very strong dependence on £. We showed its im- 
pact in the B s — > and B — > X s j decays and pre- 
sented expectations for its magnitude as dictated by RGE 
running. 

We discussed the constraints from direct searches of 
the heavy MSSM Higgs bosons. Bounds from H/A —> 
t + t~ searches mainly depend on Ma and tan/3 and 
are robust against variations of other SUSY parame- 
ters. Separately, searches in the H/A — > bb channel 
show a stronger dependence on the parameters under 
consideration, in particular on the sign and magnitude 
of the Higgsino mass parameter, /1, and therefore pro- 
vide complementary information. Currently, however, 
the H/A — > t + t~ searches are more strongly constrain- 
ing for the considered scenarios. 

On the flavor side, we considered the tree level decay 
B — >• tv as well as the loop induced FCNC processes 
B s — > /i + /i _ and B — » X s ^. The recent experimental 
updates on the BR(£? — > tv) show reasonable agreement 
with the SM prediction. At tree level, charged Higgs 
contributions to B — > tv interfere destructively with 
the SM amplitude. At the loop level, a net construc- 
tive interference is in principle possible for very large 
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and negative /xtan/3. However, we find that the cor- 
responding regions of parameter space are excluded by 
vacuum meta-stability considerations. The B — > tv de- 
cay can lead to constraints in the M^-tan /3 plane also for 
M A > 800 GeV where current direct searches for MSSM 
Higgs bosons end. For such heavy Higgs bosons however, 
B — > tv only probes very large values of tan/3 > 60. The 
B — > tv constraints depend only moderately on SUSY 
parameters other than Ma and tan fj. In particular, they 
depend only weakly on possible new sources of flavor vi- 
olation beyond the MFV ansatz. 

The constraints from the FCNC decays on the tan/3- 
Ma plane depend crucially on several parameters, in 
particular the Higgsino mass, (i, the stop trilinear cou- 
pling, A t , the gluino mass, M 3 , the mass splitting of 
the left-handed squarks, AQ13, and its alignment in fla- 
vor space, C- The current experimental bounds on the 
BR(B S — > n + ijr) lead to strong constraints in the large 
tan j3 regime of the MSSM with MFV. Constraints are 
particularly strong if the MSSM contributions interfere 
constructively with the SM, which happens for sign(fj,A t ) 
= -1. In that case, even for moderately large tan /3 <~ 30, 
heavy Higgs masses of up to 1 TeV can be probed. Note 
that these bounds can have a strong dependance on AQ 13 
and (. For negative A t , they become less constrain- 
ing for larger values of £ and larger AQ13. The main 
dependence is to a good approximation on the product 
C x AQ13. In a mSUGRA setup this product is cor- 
related with the mass splitting of the right-handed up 
squarks ( x AQ13 ~ AU13/2. From our RGE analy- 
sis of a simple mSUGRA model, we expect £ = 0.8 for 
tan/3 = 45 and ( even closer to 1 for smaller tan/3. We 
also find AQi 3 <~ 20% to 35%, which should be approx- 
imately generic for SUSY breaking models with flavor 
universal soft masses at the GUT scale and light gluinos 
M 3 < 2 TeV. Such values have visible impact on the 
bounds derived from BR(_B S — > fi + If a lower bound 
on BR(_B S — > above one half of the SM predic- 

tion is established in the future, destructively interfering 
SUSY contributions will also be highly constrained. 

It is important to stress that the MSSM contribu- 
tions to B s — > do not necessarily decouple with 
the SUSY scale, but can probe masses of SUSY parti- 
cles far above the scales that are currently reached by 
direct searches. On the other hand, the MSSM contribu- 
tions to the B — > X s -f decay do decouple with the SUSY 
scale, but even so, the B — > X s j decay can give non- 
trivial constraints on the MFV MSSM parameter space. 
If SUSY particles are heavier than <~ 2 TeV, charged 
Higgs contributions to BR(£> — > X s j) still lead to a con- 
straint for small Ma which is almost independent of all 
other parameters if tan/3 is not large. The correspond- 
ing bound in the M^-tan/3 plane can be stronger than 
the bounds from direct searches for tan/3 < 5 and rules 
out Ma < 300 GeV if squarks are decoupled. For a 
TeV scale SUSY spectrum, SUSY loops can also con- 
tribute sizably to B — > X s j. This is particularly true 
for a sizable mass splitting AQ13 and negative values of 



At, where Higgsino and gluino loop contributions add 
up constructively. Again, £ can impact the implied con- 
straints significantly. In contrast to B s — »■ how- 
ever, the bounds become stronger for larger values of (, 
if A t is negative. A main conclusion of our work is that 
the current bounds from B — > X s -y and B s — > arc 
minimized if both /z and A t are positive. In this region 
of parameter space, (g — 2)^ generically prefers a positive 
sign of M2. 

We remark that the discussed FCNC B decays are also 
sensitive to sources of flavor violation beyond MFV. For 
the MSSM with generic flavor violating structures, how- 
ever, bounds from FCNC processes become significantly 
more model dependent. 

Finally, we analyzed the impact of the updated bounds 
from dark matter direct detection searches. We found 
that the parameter space region where Mi ~ /1 is ruled 
out throughout the whole M^-tan/3 plane. Away from 
bino-higgsino degeneracy, the current XenonlOO bounds 
still give strong constraints in the M^-tan /3 plane as long 
as Mi and \x are below 1 TeV and /U is positive. The direct 
detection bounds are minimized for negative \i, where 
light and heavy scalar contributions to the neutralino- 
proton cross section partially cancel. These direct detec- 
tion constraints are the least robust among the consid- 
ered bounds, since they are subject to important nuclear 
and astrophysical uncertainties and depend crucially on 
the assumption that the lightest MSSM neutralino con- 
stitutes the entire dark matter in the universe. 

In summary, we presented the viable MSSM param- 
eter space using the MFV assumption, incorporating 
the discovery of a Higgs state at 125 GeV, the null 
direct search results for supersymmetric particles and 
for H/A — > t + t~ and bb, and constraints from B and 
K observables as well as dark matter direct detection 
searches. We also discussed and imposed electroweak 
vacuum meta-stability requirements, and we illustrated 
expectations for B flavor bounds arising from a renor- 
malization group running analysis of generic minimal su- 
pergravity models. Throughout, we have emphasized the 
connection between flavor observables and direct collider 
searches in exploring the MSSM parameter space. This 
complementarity is not only important for understanding 
the present status of the MSSM with MFV, but it is also 
central to interpreting future experimental discoveries. 
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APPENDIX: LOOP FUNCTIONS 

The loop induced "wrong" Higgs couplings involve a 
single loop function 



g{x,y,z) = 



xlogx 



+ 



y logy 



(x-y)(x-z) (y-x)(y-z) 
| z log z 
(z-x){z-y) ' 



The loop functions hj t % enter the charged Higgs con- 
tributions to the b — > sj transition 



h 7 (x) 
h$(x) 



3 — 5x 



2- 3a; 



12(1 -a;) 2 ' 6(l-x)3 l0gX ' 



3-a: 1 
4(1 - a;) 2 + 2(1 - a;) 3 



log a; . 



The loop functions that enter the Higgsino, gluino, and 
Wino contributions to the b — > s-f transition can be writ- 



ten as 

= fl + g/2 , /| = _ g/2 , ffi = —f3 - ! 
1 - 3 

fH _ f f9 — f 9f fW _ f 

.IS — /2 > J 8 — — 3/2 — •J/l , h — ~2 ' 



with 



fi(x,y,z) = - 



+ 



x 2 log x 



y 2 log y 



f 2 {x,y,z) = 



(x - y) (x - z) 3 (y-x)(y- z) 3 
(x 2 y 2 - 3a;yz 2 + (x + y)z 3 ) log z 

(x — z) 3 (y — z) 3 
x(z — 3y) + z(y + z) 
2{x - z) 2 (y - z) 2 ' 

xz log x yz log y 



+ 
+ 

h{x,y,z) = - 
+ 



(x-y)(x-z) 3 {y-x)(y-z) 3 
z(xy(x + y) — 3xyz + z 3 ) log z 

[x — z) 3 (y — z) 3 
z(y - 3z) + x(y + z) 



2(x-z) 2 (y-z) 2 ' 

z 2 log x z 2 log y 



(x-y)(x-z) 3 (y-x)(y-z) 3 
z 2 (x 2 + xy + y 2 - 3(x + y)z + 3z 2 ) \ogz 

(x — z) 3 (y — z) 3 
x(y — 3z) + z(5z — 3y) 
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